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Abstract. We consider the problem of determining the minimal time for which an energy 
supply source should operate in order to supply a system with a desired amount of energy 
in finite time. 



1. Introduction 



While boiling water or any other liquid, most of us have noticed that the heater can be 
switched off at an intuitively chosen time and the liquid will still reach its boiling point no 
long after the heater is switched off. A natural question arises: how can switch-off time be 
chosen in an optimal way so that electrical energy will not be wasted? In other words what 
is the earliest time at which the heater can be turned off while still reaching the liquid's 
boiling point in a finite time? A general formulation of the problem is as follows: Let D be 
a device that supplies energy to a system S through a supply line. Let E'(t) be the energy 
supply rate. We assume that D can be switched on and off and that it continues to supply 
energy, at a decreasing rate, for some time after it has been switched off. Question: What 
is the minimum switch-off time of D (corresponding to the minimal operational time of D) 
in order to transfer to S a total amount of energy Q (where Q > is given)? The device D 
can also be viewed as a control mechanism for bringing the system S from an energy level 
El to a higher energy level E2 in finite time while operating for the minimum time possible. 
Clearly, the solution of this problem can have a lot of applications , both civilian (e.g energy 
conservation) and military (e.g minimizing detection risk). 

2. Examples 

EXAMPLE 1. Exponential Model. 

We consider a simple example of an energy device supplying energy to a system at time t > 
at a rate (which in what follows we consider to include the rate at which the transferred 
energy is possibly radiating from the system and/or the supply line) given by 



where a and b are positive real numbers characteristic of the source but also depending on 
the environment, to is the time at which the energy supply rate is at its peak, and ti is 
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the switch-off time. We assume that the optimal switch-off time is after the rate of energy 
supply has been stabihzed i.e that ti > to. We also assume that after switching-off at time 
^1 ^ ^0) the source stops transferring energy to the device at time ti + T, where T > is 
independent of ti. Let Q > Ohe the amount of energy that we wish to transfer to the system. 
We assume that the transfer of this amount of energy will occur at some time t2 > h > tg. 
The energy Q could be, for example, the energy required to bring a liquid substance to its 
boiling temperature, or the energy required for complete phase transition. In the latter case 
Q = rriLy, where m is the mass of the liquid and Ly is its latent heat of vaporization (cf. 
[2]). We require that 



t2 



E'{s) ds = Q 



which implies that 



/ E'{s)ds+ I E'{s)ds+ I E\s)ds = Q 

Jo Jto Jti 

or, by the definition of 

/ (e- -l)ds+ (e"*° -l)ds+ J ^ (e-^(-*i) - e''^) ds ^ Q 

Jo Jto Jti ^ 



which implies that 



— s 



s=to 



s=0 



+ (e«*°-l) {h-to) + f 



^ato 



>(s-ti) 



S=t2 



S=tl 



Q 



from which letting 

y 

Lo 



L2 :-- 



t2 — ti 

1/1 p«*o 1 _ pato 1 

e«*o -1 \ a a 1 - e"^^ b 

1 



6(1 - e-b^) 
1 



1 - e- 



-bT 



we obtain 



(2.1) 



Notice that if |/ = then = Lq -|- Li is the time the device supplies the desired energy level 
Q without being switched-off. The maximum value of y is y = Umax — {h + T) — ti — T . 



MINIMAL OPERATION TIME OF ENERGY DEVICES 3 

The optimal switch-off time ti is therefore determined from (12.11) by letting y = T and it is 
given by 



fi = Lo + Li e-^^ + L2 T. 
Energy level Q is reached at time 

t2 = ti + T = Lo + Li e"^^ + (L2 + 1) T. 
EXAMPLE 2. Linear Model. 

A simplified version of the previous model is obtained by assuming that 



— f 

to ^ 



no <t<to 

E'it) = { a ifto<t< ti 

-f(t-ti)+a ift>ti 

where to > 0, a > 0, and T is as in Example 1 above. In this case 



t2 



E'{s) ds = Q 



implies that 



to 



E'{s)ds+ I E'{s)ds+ I E\s)ds = Q 
or, by the definition of E'{s), 



to 



t2 



I —sds+ ads+ i——{t — ti) + a] ds = Q 
Jo io Jto Jti ^ T / 



which implies that 



a s 



S=to 



+ a{ti - tQ 



'-0 -Js=0 

and letting y = t2 — h we obtain 



-—{s-ti) +as 



-t2 



Q 



. 1 2 , Q to 
ti = 777fy -y + - + iz- 

21 a 2 

As in Example 1, substituting y by ymax = T we obtain 

1 ^ Q ^0 
ti = — T + - + - 

2 a 2 

which is bigger or equal to to if and only if^ — |T — y>0. Energy level Q is reached at 
time 
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^ _ to Q 1^ 
t2 = ti + T = ^ + ^ + -r. 

^ (X ^ 

3. General Description of the Optimal Switch-Off Time 
The examples treated in detail in the previous section suggest the following general theorems. 

Theorem 1. Let, for each ti > to 

( fit) ifO<t<to 
Kit) = < fito) ifto<t<t, 
[ 9{t) ift>t^ 

where f is continuous and increasing with f{0) = 0, and g is continuous and decreasing with 
gifi) = f{to)- Let F and G denote the anti- derivatives of f and g respectively. If there exists 
T > such that g{ti + T) — for all ti >to, then the optimal switch-off time ti > to is the 
solution of 

F{to) + fito) it\ - to) + G{t\ + T) - G{t\) = Q. 

Proof. The condition 

r E',^{s)ds^Q 
Jo 

implies that 

{F{to) - F{0)) + fito) {ti - to) + {G{t2) - G{ti)) = Q 
which by F(0) =0 and the fact that for the optimal ti, t2 = ti + T implies that 



F{to) + fito) iti - to) + Giti + T) - Giti) = Q. 



□ 



We can generalize the above theorem to the case of a heat source that supplies energy at a 
possibly strictly increasing (or even arbitrary continuous) rate as follows. 

Theorem 2. Let, for each ti > to, the energy function E^^it) be an increasing continuously 
differentiable function of t. The optimal switch-off time ti and the associated time t2 at 
which the energy level reaches Q are determined by the system of equations 

Ef^it2)^Q, E'M^Q. 
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Proof. It is clear that the optimal switch-off time ti makes full use of the energy source in 
the sense that the desired energy amount Q is supplied at the moment when the source dies 
out i.e when E'- (^2) = 0. For such a ^2, 

^1 




E'^-^{s) ds = Q 



implies 

f ' E[^^ (s) ds+ r E[~^ (s) ds = Q 
Jo ^ Jt\ ^ 

i.e E^^{ii) - E^^{0) + ^4(^2) - ^4(^1) = Q and since ^<;(0) = 0, we obtain E^^{t2) = Q. 

□ 

4. Noisy Supply Line 

The energy supply rate functions := E'{t) used in Examples 1 and 2 above, can be 

viewed on each of the intervals [0,^0]; [^Oi^i] and [^1,^2], as solutions of ordinary differential 
equations of the form 



(4.1) d^tAt) = ici(ptAt) + C2)dt 

for appropriate constants ci,C2 (depending on the interval) and initial condition 0ti(O) = 0, 
with the different branches of (f)ti{t) tied up at to and ti. It could be the case however, 
that the supply line connecting the energy source D with the system 5" is affected by noise, 
appearing in (14.1 p in the form of additive noise 



(4.2) d(t)t, (t) = (ci (f^t, (t) + C2) dt + (C3 (Pt, (t) + C4) dB{t) 
or equivalently 

(4.3) M)= I {ci(f)tM + C2) ds+ [ {C3(l)t,{s) + C4) dB{s) 

Jo Jo 

where Ci, C2, C3, C4 are constants, B{s) is one dimensional Brownian motion and the stochastic 
integral on the right hand side of (14.31) is in the sense of Ito (cf. [1]). In that case 4>ti(t) is 
actually a stochastic process (pn {t, i^) and the problem of finding the optimal switch-off time 
ti now amounts to finding the first ti for which 




fitiis) ds = Q 



for some finite t2 > ti, where ^tx{t) denotes the mathematical expectation of 0ti(t, uj) . Since 
the mathematical expectation of an Ito stochastic integral with respect to Brownian motion 
is equal to zero, (14.31) implies upon taking the expectation of both sides that 
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/"ii (t) = / (ci /Ufi (s) + C2) (is 
which can be solved exphcitly and yields the mean energy supply rate 

Cl 

where on each interval the constants Ci and C2 are determined by using the initial and tying 
up conditions. Thus we are reduced to the deterministic problem considered in the examples 
of Section 2, but this time for the mean energy supply function. The method extends directly 
to the case when equations (14.11) and (14. 2 p are replaced by the more general equations 

and 

d(pt, (t) = fit, (pt, (t)) dt + g{t, <Pt, it)) dB{t) 

respectively. 
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